Ferromagnetic Kondo effect in a triple quantum dot system 
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We propose that a simple device of three laterally-coupled quantum dots, the central one contacted 
by metal leads, can realize the ferromagnetic Kondo model, which is characterized by interesting 
properties like a non-analytic inverted zero-bias anomaly and an extreme sensitivity to a magnetic 
field. Furthermore, by tuning the gate voltages of the lateral dots, this device may allow to study 
the transition from ferromagnetic to antiferromagnetic Kondo effect, a simple case of a Berezinskii- 
Kosterlitz-Thouless transition. We model the device by three coupled Anderson impurities that 
we study by numerical renormalization group. We calculate the single-particle spectral function of 
the central dot, which at zero frequency is proportional to the zero-bias conductance, across the 
transition, both in the absence and in the presence of a magnetic field. 



PACS numbers: 75.20.Hr, 73.63.Kv, 73. 21. La 

In spite of its simplicity - a magnetic impurity embed- 
ded in a conduction bath - the Kondo model exhibits rich 
many-body physics [1 that continues to attract scientific 
interest in a variety of contexts. It is common to dis- 
tinguish between Kondo models with Fermi-liquid prop- 
erties a' la NozieresL2 from those that instead display 
non-analytic, hence non-Fermi-liquid, behavior as func- 
tion of state variables, and which include under- and over- 
screened Kondo models [2 [3111] as well as clusters of mag- 
netic impurities in particular circumstances [5 - 7J. Non- 
Fermi liquid properties are not common in traditional 
magnetic alloys, [8] where the metal hosts generally pos- 
sess enough scattering channels that can perfectly screen 
the magnetic impurity. They may instead be realized in 
confined scattering geometries such as a quantum dot or 
a magnetic atom/molecule contacted by metal leads. In- 
deed, by means of such devices there are already many ex- 
perimental realizations of exotic non-Fermi-liquid Kondo 
models, see e.g. Refs. 0, [TO], [TT] and [H]. 

One case, however, which so far remains elusive is the 
ferromagnetic Kondo model (FKM)[13j - where the im- 
purity and the conduction electrons are coupled ferro- 
magnetically - except for its indirect manifestation in 
the under-screened Kondo effect. O [1] It has been pro- 
posed that the so-called giant moments induced by 3c? 
transition metal impurities diluted in 4c? transition met- 
als may actually be a manifestation of FKM,[T3] but ex- 
periments that could pin it down are still lacking. This 
is unfortunate since the FKM is the simplest example of 
non-Fermi liquid behavior; at low temperature the impu- 
rity spin behaves essentially as a free local moment apart 
from logarithmic singularities. [TOl [T5] 

Here we present a possible realization of a FKM by 
means of three laterally-coupled quantum dots. We 
also discuss the appealing possibility of crossing the 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition 
from ferromagnetic to antiferromagnetic Kondo effect, 




FIG. 1. (Color online) A schematic representation of our de- 
vice described by Eq. [ij with three quantum dots (in green). 
Only the central one, labelled as 0, is attached to metallic 
leads (in red). 

whose spectral weight anomaly change we study here by 
means of numerical renormalization group (NRG) [lBl418| 
in a toy-model for the device. 

Our gedanken (but entirely feasible) T9l [20] set-up, 
schematically shown in fig. [T] consists of three quan- 
tum dots, labelled as ±1 and in the figure, with the 
central one contacted by two metal leads, R and L. We 
model the isolated three-dot device with a Hamiltonian 




- E i^- ^-1-^0. + U cLc+i + H.c.) , (1) 

where we keep just one orbital per dot, Ui are the 
charging energies, and the dots are mutually coupled by 
single-particle tunneling. As usual, we shall assume non- 
interacting leads, coupled to the central dot by tunnel- 
ing. For convenience, we also assume equivalent R and L 
leads , so that only their symmetric combination matters 
in the linear response regime of interest to us. Particle- 
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hole symmetry is also assmxied. The lead-dot tunneling is 
thus parametrized by a single quantity, the hybridization 
width r = TT^jj. |Vk|^ S{e]^), where is the energy and 
Vk the tunneling amplitude into the dot of the symmetric 
L + R combination of the lead electrons at momentum k. 

We further assume that each dot is brought by gate 
voltage in the Coulomb blockade regime with a single 
unpaired electron, i.e. ~ and t± <C Ui, \fi. In this 
limit, the isolated trimer behaves like a three-site Heisen- 
berg model described by an effective Hamiltonian 



U 



J+ So • S+i + J_ So • S_i, 



(2) 



with positive J±, where S^, i = 0, ±1, are spin-1/2 op- 
erators residing on the corresponding dots. The ground 
state of ([2]) has total spin 1 /2 and explicitly reads 



GS, cr) =cos6i I O, cr) - si 



E,a), 



(3) 



where a J, is the z-component of the total spin, and 
tan26l = v^(j+ - J-)/(J+ + J-). In Eq. (B, | 0,ct) 
is the state, odd by inversion through dot 0, obtained by 
coupling dots -f 1 and —1 into a triplet and coupling the 
latter to dot to form a spin 1/2. Vice versa, | E,(t) is 
the state, even by reflection, obtained by coupling dots 
-1-1 and —1 into a singlet, leaving a free spin-1/2 on dot 
0. 

If on dot the electron is removed or one more elec- 
tron is added through the leads, the trimer ends up in a 
triplet or singlet configuration, with probability propor- 
tional to cos^ 9 and sin^ 9, respectively. Specifically, if 

^ ~ y^X^k l^kl^ ^ Uq, the Kondo exchange, Jeg, can be 
found by second order perturbation theory: 
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and where 
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are the Sz = component of the triplet state and the 
singlet state, respectively. It follows that 

cos^ 9 
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where the resolvent operator 
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and 7s > 7t > 0, because the intermediate singlet state 
has lower energy than the triplet. The lead-dot exchange 
is therefore ferromagnetic if 7* cos^ 6* > 3js sin^ 6*, and 
antiferromagnetic in the opposite case, while 74 cos^ 9 — 
3 7s sin^ 9 actually identifies the transition between the 
two. We observe that, if inversion symmetry holds, 
J+ = J_, then 9 — hence the lead-dot exchange is 
ferromagnetic, providing therefore a realization of the 
FKM. We expect that in a real device inversion sym- 
metry will be generally broken; nevertheless there still is 
a good chance for ferromagnetism to survive in a wide 
region (by definition cos^ 9 > 3 sin^ 9, hence just because 
7s > 7t is it possible for the Kondo exchange to turn 
antiferromagnetic) . 

In conclusion, the set-up shown in fig. [T] seems indeed 
able to realize the much-sought FKM. Moreover, it sug- 
gests a simple way to study experimentally the transition 
from the FKM to the more conventional antiferromag- 
netic Kondo model, first described by Anderson, Yuval 
and Hamann [5T], and expected to be of the BKT type. 
Indeed, changing the gate voltage e+i with respect to 
e_i drives the system further away from the inversion 
symmetric point, eventually turning the exchange from 
ferromagnetic to antiferromagnetic. We investigate this 
possibility by considering the Hamiltonian ([l]) in the sim- 
ple case when t+ = t- = t, Uq = U+i = C/_i = U, 
and studying it as function of e+i = — e_i — Se by 
means of NRG. fTBHTS] For simplicity we take a flat 
conduction-band density of states, p(e) = po — 1/2D 
when e e [—D,D] and zero otherwise, with the half- 
bandwidth D our unit of energy. Calculations employ 
a home-built NRG code, which implements as quantum 
numbers charge and the z-axis projection of the total 
spin. We use the value A = 2 for the discretization pa- 
rameter of the conduction band, and keep 300 states per 
subspace, which roughly means 3500 total states per it- 
eration. The spectral functions (i — —1, 0, +1) 
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are computed through the patching technique, 
and delta-peaks are broadened in the form 



with b — 0.7. The temperature is set by the length of the 
NRG chain and is almost zero, T ~ 10"^!?. 

In fig. [2] we show the single-particle spectral density 
on the central dot, Aq{uj), as a function of the frequency 
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FIG. 2. (Color online) NRG spectral function Ao{uj) on dot 
as a function of Se for U = 0.3, t = 0.03, T = 0.02 tt with 
a flat conduction band density of states of half-bandwidth 
D = 1. Note the transition from the ferromagnetic Kondo, 
signaled by a minimum of Ao{uj) at = 0, to the regular, 
antiferromagnetic Kondo, where instead Ao{(jj) is maximum 
at a; = . 

Lu for different values of Se. Aq{ll! — 0) is actually propor- 
tional to the zero-bias conductance. j24] while its behav- 
ior at small uj mimics that of the conductance at small 
bias. We note that for small 6e, which corresponds to 
a ferromagnetic Kondo exchange, Jeff < in Eq. 
Ao{u} = 0) = 0; the zero-bias conductance vanishes and 
the zero-bias anomaly is inverted. This spectral func- 
tion displays a logarithmic "dimple" at w = 0, asso- 
ciated with the 1/ln^ (w/Tg) singularity typical of the 

FKM.dS] Here Tq oc ypoffi exp {-l/po Jeff) is an en- 
ergy scale that actually diverges on approaching the BKT 
transition from the ferromagnetic side, i.e. J^s — > from 
below. [13^ 

On the contrary, above a critical Sec, Aq{uj = 0) sud- 
denly jumps to a finite value; the conventional Abrikosov- 
Suhl spectral anomaly of the regular, antiferromagnetic 
Kondo model is recovered. Because of our choice of 
parameters, the Hamiltonian is particle-hole symmet- 
ric hence Ao(a; = 0) jumps from to its unitary value 
2/(7rr) at the transition. Away from particle-hole sym- 
metry, the jump still exists but will be smaller. The 
energy scale that controls the antiferromagnetic side, 
Jeff' > 0, is of course the Kondo temperature Tk oc 
exp ( — 1 /pq Jeff) , usually defined as the half-width 
at half- maximum of Ao(w), which vanishes on approach- 
ing the BKT transition. Jog' — > from above. 

The standard, most reliable way to reveal the Kondo- 
like origin of a zero bias anomaly is by applying a mag- 
netic field B. In the standard antiferromagnetic Kondo- 
effect, a magnetic field will split the Abrikosov-Suhl res- 
onance only if sufficiently large, gfisB > 0.5 T^f pS]. 
This is indeed the case on the antiferromagnetic side of 




FIG. 3. (Color online) NRG spectral function Ao{uj) of the 
central dot 0, for increasing magnetic field B (we take fiB = 1 
and g — 2, hence 2B is the Zeeman splitting) and the same 
parameters of Fig. |2] Panel (a) refers to Se = 0, while 
panel (b) to Se = 0.0147. In the ferromagnetic Kondo (FMK) 
regime, panel (a), the logarithmic dimple is immediately de- 
stroyed, and replaced by inelastic spin-flip excitations; in the 
antiferromagnetic (standard) Kondo regime, panel (b), where 
Tk ~ 6 X 10"^ , the Abrikosov-Suhl resonance is only split 
by a sufficiently large field 2g^BB ~ ksTK ■ 



the transition, Jeff > 0, see panel (b) of fig. |3] On 
the contrary, on the ferromagnetic side, Jeff < 0, panel 
(a), any magnetic field, however small, destroys the loga- 
rithmic dimple replacing it right away with a symmetric 
pair of inelastic spin-flip Zeeman excitations. In addition, 
Ao{uj = 0), hence the zero-bias conductance, increases 
with B at low temperature, contrary to the antiferro- 
magnetic side, where it drops. We expect moreover that 
a finite temperature T will cutoff the logarithmic dimple 
at low frequency and raise up Ao{uj = 0) ^ l/ln^(r/To), 
thus leading to an increase of zero-bias conductance, 
again unlike the regular antiferromagnetic Kondo effect. 

In conclusion, we have shown that a three dot device, 
the central one contacted by metal leads, may provide 
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a realization of the ferromagnetic Kondo model. Gating 
of the lateral dots will in addition drive a Berezinskii- 
Kosterlitz-Thouless transition from ferromagnetic Kondo 
to regular, antiferromagnetic Kondo effect. The two 
phases should differ sharply in their zero-bias conduc- 
tance anomaly, the ferromagnetic one being inverted and 
very differently modified by a magnetic field and by tem- 
perature. 

More generally, our proposed system illustrates a 
generic mechanism leading to FKM, namely tunneling 
across one orbital in presence of other magnetic or- 
bitals.This kind of situation could for example also be 
realized at selected surface adsorbed molecular radicals 
and detected in, e.g., STS or photoemission anomalies, 
an area where there is much active work. However, we 
should mention that, according to our calculations, the 
FKM anomalies are more visible when both U and the 
tunneling amplitude into the leads are larger than the 
inter-dot tunneling t, which might be hard to achieve in 
molecular radicals. 
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